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Abstract: A binary operator * over real numbers is said to be associative if 
(x * y) * z = x * (y * z) and is said to be reducible \ix*y = x*z or y*w = z*w 
if and only if z = y, The operation is said to have an identity element e if x * e = x. 
In this p aper a characteriza t ion of a subclass of the reversed generalized Pareto dis- 
tribution I Castillo and Hadi 1 1995)) in terms of the reversed lack of memory property 
• lAsha and R ejeesh (2001])) is generalized using this binary operation and probability 

distributions are characterized using the same. This idea is further generalized to the 
bivariate case. 
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O . 1. Introduction 



A binary operation * over real numbers is said to be associative if 



(x * y) * z = x * (y * z) (1-1) 

a | 

for all real numbers x,y, z. The binary operation * is said to be reducible ifx*y = x*z 
if and only if y = z and if y * w = z * w if and only if y = z. The general reducible 
continuous solution of the functional equation ( 11.11 ) is ( Aczel , 19661 p. 254) 



x*y = g 1 {g(x)+g(y)) (1.2) 

where g(-) is a continuous and strictly monotone function provided x,y,x * y belong 
to a fixed interval A in the real time. The function g(-) in dl.2t is determined up to a 
multiplicative constant; for all x, y in a fixed interval A, implies g^ = agi for all x 
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in that interval, for a ^ 0. We assume hereafter that the binary operation is reducible 
and associative with the function g(-) continuous and strictly increasing. Further more 
assume that there exists an identity element e 6 R such that 



x * e 



x G A. 



(1-3) 



Further more every continuous, reducibl e and assoc iative operation defined on an inter- 
val A in the real line is communic ati ve dXczell d 1 966l)) . Character izati on of distributions 

throug h bi nary operations is gi ven in lMul iere and Scarcini ( 1987 ) and Mul iere and Prakasa Raol 
(2003). In Prakasa Rao (2004), the bivariate lack of memory property (Roy (2002)) is 
generalized and classes of bivariate probability distributions which include bivariate 
exponential, Weibull, Pareto distributions are characterized under binary associative 
op erations. 

lAsha and Reieeshj ( 20071) character ized a subclass of the reversed generalized Pareto 
distribution (Castill o and Hadil ( 1995)) using the reversed lack of memory property. In 
this paper we generalize the reversed lack of memory property and characterize prob- 
ability distributions using this property. In Section 2 the univariate reversed lack of 
memory property is generalized using the binary operation and a class of probability 
distributions which include subclass of reversed generalized Pareto distribution, power 
function and reflected Weibull is characterized. In Section 3 this property is extended 
to the bivariate set up the class of bivariate probability distributions which are charac- 
terized by this property is derived. 



2. Univariate characterizations 



Let X be an absolutely continuous random variable with density, distribution function 
and survival function defined as /(•), F(-) and R(-) respectively. Let a — mi{x : 
F(x) > 0} and b — sup{x : F(x) < 1}, then (a, b), — oo < a < b < oo is the interval 
of support of X. A random variable X with support (a, b) where a < and b > is 
said to have the reversed lack of memory property (RLMP) (Ash a and Rejeeshl d2007l) 1. 
if 

P(X < x\X < x + t) = P(X < 0\X < t) 

for all x and t such that a<x<x + t<b. 

In terms of the distribution function we can write the RLMP as 



where a < x < x 



F(x) ■ F(t) 
t < b with a < 0. 



F(x + t)-F(0) 



(2.1) 



Solving (El) (L\czelLll966[ p. 38) we have 



F(x) = 



c(x — b) 



< b 

> 6, c> 0. 



(2.2) 



This type 3 extreme value distribution also belongs the the reversed generalized Pareto 
distribution developed by lCastillo and Hadil (11995b as a fatigue model that satisfy cer- 
tain compatibility conditions arising out of physical and statistical conditions in fatigue 
studies. 
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We shall extend the RLMP using the binary operator * as 

P(X < x\X < x*t) = P(X < e\X < t) (2.3) 
for all a < x < b, a < e, e < t < b, x * t < b. 

Here e is specified as in ( 11.3l l. The general reducible continuous solution is (11.2b . 
In terms of the distribution function ( 12.31 can be written as 

F(x)-F(i) = F(x*i)-F(e) (2.4) 

for all a < x < b, a < e, e < t < b, x * t < b. 

We refer to ( 12.4b as generalized reversed lack of memory property (GRLMP). 

Recently the reversed hazard function ha s become quite popular since it is very 



useful in the analysis of left censored data (see. lGupta and Hannl(l200ll) . lAndersen et al 



(119931) '). The reversed hazard rate (RHR) function of X is defined for x > a as 



It has been observed that there exists a relation between proportional reversed hazard 
class of distributions and the exponentiated class of distributions. The class of dis- 
tributions [F(-)] Q , a > can be defined as the exponentiated class of distributions 
with base distributions F(-). In this case the reversed hazard rate function of X satisfy 
\ Y (t) = a\ x {t) where Y has a distribution [F(-)] a . Observe that F(x) in d2~2l is a 
proportional reversed hazards model. 

In the next theorem we show that the continuous solution of d2.41 > are generalized 
proportional reversed hazard (PRH) models. 



Theorem 2.1. The (continuous) solution of (12.41 ) is 

F(x) = exp[c(<7(a;) — g(b))\ with c > and x 6 (g~ 1 (a),b). 
where g is a continuous and strictly monotone function. 
Proof. Substituting i ll .21 in ( 12.41 ) we have 

F(x) ■ F(t) = F[g-\g(x) + g(t))} • F(e). (2.5) 



s = g(x), u = g(t) and F o g^ 1 = H (2.6) 



Writing 
( 12.51 ) becomes 

H{s) ■ H(u) = H{s + u) ■ H(g(e)) 
for all g(a) < s < g(b) and 17(e) < u < g(b) which implies 

H(s) 



G(s) ■ G(u) = G(s + u) where G(s) 



H{g{e)) 



for all g(a) < s < g(b) and g(e) < u < g(b) which is the Cauchy functional equation 
and the solution is 



G{s) = e cs , s G (g(a), g(b)) and c> (lAczellll966l p. 38). 
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From ( I2.6l l we get 

F(x) = e c ^ ■ F(e), x £ (g-\a),b). (2.7) 

Now taking x = b we get F(e) — e~ cg ( b \ 
Thus (12.7b reduces to 



F(x) 



' e c[g(x)-g(6)] ) ^ g ( 5 ~1( ) ) fc) 

1, x > 6, c > 0. 



Hence the theorem. □ 
If we particularize the operation *, we obtain different types of di stributions. 



Example 1. Subclass of reversed generalized Pareto distribution (ICastillo and Hadi 



dl995h l. 



For x * y — x + y, we get g{x) — x, x € (— oo, b) with a — — oo, b < oo and e = 0. 
The distribution function is now given by 



F(x) 



exp[c(x — b)], x < b 

1, x > 6, c > 



Example 2. Power function distribution. 

(i) For x * y = xy, we get g(x) = logic, x G (0, 6) with e = 1. In this case the 
distribution function is given by 



F(x) 



(f) c , 0<x<b 
1, x > 6, c > 0. 



which is the distribution function for the power function distribution. 

(ii) If we take x*y = x + y + xy, we get g(x) = log(x + 1), x S (—1, 6), 6 < cxd 
with e = 0. Here the distribution function is given by 

F(x) = lmr> -i<*<b 

w 1 1, x > b, c> 0. 

which is the distribution function for the power function distribution in the support of 

(-l,6),6<oo. 

Example 3. Reflected Weibull distribution dLai and Xie 



If we take x * y = \J x 1 + y 2 , we get g(x) — —x 2 , x € (— oo, 0), with e = 0. The 
corresponding distribution function is given by 



e~ cx , x < 
1, x > 0, c> 



F(x) 

which is the distribution function for a reflected Weibull distribution. 
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Remark 2.1. If X is a random variable in the support of (g (e), b) with F{e) ^ 0, 
then X satisfies the GRLMP if and only if 



F x (x) 



exp[c[ 5 (a;) - g(b)]], g 1 (e) <x<b 
1, x > b, c> 



which has a probability mass at x — g 1 (e). 

Now we try to evolve the above concept to the higher dimensions. One of the main 
problems associated with such an attempt is that there is no unique way of evolu- 
tion. We here consider th e definition of bivariate reversed lack of memory property 
( Asha and Rejeeshl (2007)) to extend the concept of univariate GRLMP to the bivari- 



ate case and derive bivariate models characterized by the respective property. Since 
multivariate derivations are rather straight forward extensions it is excluded. 

3. Bivariate Extensions 

Consider a random vector X = (Xi,X2) in the two-dimensional space with joint 
distribution function F(xi,x 2 ) = P(X\ < x%,X 2 < x 2 ) in the support of (o, b) 2 
where 

a = mf{xi\Fi(xi) > 0} and b — sup{xi\Fi(xi) < 1}, 

with F(xi, b) and F(b, X2) as the marginals of X*, i = 1,2. 
Then a direct extension of RLMP is 

F(x 1 +t 1 ,x 2 +t 2 )-F(0,0) =F(x 1 ,x 2 )-F(t 1 ,t 2 ) (3.1) 

for all Xi and t{ such that a < Xj < Xj + U < b, a < 0, i = 1, 2. 
The only solution for ( 13. Il l is ( Aczel ( 19661) ) 



F( Xl ,x 2 ) = e M*i-b)+A 2 (x 2 -b) )a< x . < 5j Aj; > 0, i = 1,2. 
We consider analogous equation of (13.1b for the associative binary operator * given by 

F(xi *h,x 2 * t a ) ' F(e, e) = F(x u x 2 ) ■ F{t x ,t 2 ) (3.2) 

for aH a < Xi < b, e < ti < b, a < e, Xi * ti < b, i = 1,2, 
Combining ( 11.2b and (13.2b we have 

F[g- 1 {g{x 1 )+g{t 1 )),g~ 1 (g{x 2 )+g(t 2 ))]-F{e,e) = F{x 1 ,x 2 )-F(t 1 ,t 2 ). (3.3) 

Writing s t = g(xi), m = g(U), i = 1,2 and F(g~\-), g~\-)) = H(-,-), 
( 13.3b becomes 

H(si + Mi, s 2 + u 2 ) ■ H(g(e), g(e)) = H(si, s 2 ) ■ H(m,u 2 ) 

for all g(a) < s, < g(b) and g(e) < Ui < g(b), i = 1,2, 
Thus 

G(si + u\, s 2 +u 2 ) = G(si, s 2 ) ■ G(ui,u 2 ) where G(si,s 2 ) - ' *' 



H(g(e),g(e)) 
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for all g(a) < s, < g(b) and g(e) < < g(b), i = 1,2. 

The solution to the above Cauchy functional equation is dAczel 



G(s 1 ,s 2 )=e x ^+ X2S2 ,X 1 ,X 2 >0 

which implies 

for all xi e (ff -1 (a),&), A, > 0, i = 1,2. 

Thus, ATi and X 2 are independent with marginal distribution functions specified by 

F Xl (xi) = e^W-O-ffW] and F X2 (x 2 ) = e ^[9(^)- 9 (b)] 

We now consider a meaningful way of extending the RLMP specified by ( 12.1b to the 
bivariate case as 

F(xi + t,x 2 +t)- F(0, 0) = F(n , ar 2 ) • F(t, t) (3.4) 
for all a < < cc^ + t < b, i — 1, 2 with a < 0. 



Then the unique solution of ( 13.4-b among distribution function is dAsha and Reieesh 



(2007)) 



F{x\ X2) = e Xl ^ Xl ^ b ^ +X2 ^ X2 ^ b ^ +Xl2 ma:x '^ Xl ^ b ' X2 ^ b ^ 
for all a < Xi < b; A,, Ai 2 > 0; i = 1, 2. 

Generalizing equation ( 13.4b using the operation *. We define the generalized bivariate 
reversed lack of memory property (GBRLMP) as 

F( Xl *t,x 2 *t)- F(e, e) = F( Xl , x 2 ) ■ F(t, t) (3.5) 

for all a < Xi < b, e < t < b, Xi * t < b, a < e, i = 1,2. 

Theorem 3.1. Let F(x\ * t, x 2 * t) ■ F(e, e) = F{x\,X2) ■ F{t, t) for all a < Xi < b, 

e < t < b, x i * t i < b, a < e; i = 1 , 2 and 

F t (x*t)-F(e) = F(x)-F l (t) (3.6) 

for all a < x < b, e < t < b, x * t < b, a < e; i — 1,2. Then the continuous solution 
of the equations (13.5b and ( 13.6b is 

F(x 1 ,x 2 ) = exp{Ai[g(a;i) - g(b)} + X 2 [g{x 2 ) -.9(6)] 

+ Ai 2 max[g(a;i) - g(b),g(x 2 ) - g(b)], 

g^ 1 (a) < Xi < b; X, > 0; i = 1, 2, A i2 > 0. 
Proof. In ( 13.5b let 2:1 = x 2 = x, then 

F(x * t, x * t) ■ F(e, e) = F(x, x) ■ F(t, t). (3.7) 

Combining ( |1.2b and (13. 7b we get 

Flg-^gix) + g(t)), g-^x) + g(t))} ■ F(e, e) = F(x, x) ■ F(t, t). 
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Taking g(x) — s, g(t) — u and F{g~ 1 {-),g~ l {-)) = H(-, •). Thus ( I3.7l i reduces to the 
Cauchy functional equation 



G(s + u) = G(s) ■ G{u) (3.8) 

H(s,s) 
H(g(e),g(e))- 



for all g (a) < s < g(b) and g( e) < u < g(b) where G(s) — — g ( s ^) 



Solving for (EH) d AczelL 1 1 966L p. 38) we have 

G(s) = e ks for k > 0. 
Thus we get pfex = e kg ^ x \ k > which implies 

F(x,x) = e k[9(x) ~ 9(b) \ x G (ff-^a)^),* > 0. (3.9) 

Hence 

F(xi,b) ■ F(t,t) 



F(xi *t,b*t) 



F(e,e) 
Fi(xi) • F{t,t) 



F(e,e) 

From ( I3.91 l we get 

e A 1 [g( a;i )- g (b)]+fc[ g (t)- 9 (6)] 



e -fcg(b) 

— e Ai[s(xi)-a(6)]+fcff(t)^ (3.10) 



Now let Xi * t = s and 6 * t = u, then by dl.2l ) we get <?(a:i) = g(s) — and 
5(6) = g(u) — g(t). Thus (13.101 1 becomes 

F( S , U) = e Al l9(s)-g(u)]+k[g(u)-g(b)] . u > g 

Arguing similarly we can prove that 

F(S, U) = e ^[9(«)-9(s)]+fc[9( S )- 9 (6)] . u < g 

Writing A12 = fc — Ai — A2 and rearranging we obtain, 

F(s,u) = exp{A 1 [ ff ( S ) - 3(6)] + X 2 [g(u) - g{b)\ 

+ A12 max[g(s) - ,g(6), ,g(u) - g(b)] 

or 

^(^1-^2) = e Al ' 5 ^ 1 - l_5 ^ b '] +A2 [ 5< ' :!:2 - l_5< ' b '] +Al2 max l9( x i)-g( b ),g{ x 2)-g(b)] ^ 

g _1 (a) < < 6; A, > 0; i = l,2andAi 2 > 0. (3.11) 

□ 
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Remark 3.1. It can be seen that a necessary and sufficient condition for X\ and X 2 
to be independent is \\ 2 = 0. 

Remark 3.2. If (X\, X2) is distributed as in Theorem \3.1\ then the distribution func- 
tion o/max(Xi, X2) has the following form: 

P[m&x{X 1 ,X 2 ) <s}= P[X X < s,X 2 < s] 
= F(s,s) 

= e (Ai+A 2 +A 12 )(9( S )-g(6))^ 

Hence, the distribution o/max(Xi, X 2 ) has the same form (with different parameter) 
as the marginal distribution of X\ and X 2 . 

Theorem 3.2. Let (Xi, X2) be a bivariate random variable with joint distribution 
function F{x\, x 2 ). Then F is distributed as in ( 13.1 11 1 if there exist an independent and 
identically distributed random variables U, V, W whose marginal distributions are 
given by F[x) = e ^3( x )-g(b)] such that Xl = m ax(Z7, W), X 2 = max(F, W). 

Proof. 

F(x u x 2 ) = P{X X < X!,X 2 < x 2 ) 

= P(U < xx, W < Xl , V <x 2 ,W < x 2 ) 

= P(U < xi) ■ P(V < x 2 ) ■ P(W < max(xi, x 2 )) 

— e Ai[90i)-g(&)]+A 2 [g(iE2)-sO)]+Ai2 max[j(ii)-j((i), g(x 2 )-g(b)} 

□ 

The distribution function F(x\,x 2 ) specified in Theorem 13. II is not an absolutely 
continuous distribution in the support of (g~ 1 (e), b) 2 with F(e, e) 7^ 0. In this case we 
can write 

F{x l ,x 2 )=F d + F c + F ac 



where 



e Ai(e-g(f))). 



g 1 (e), x 3 ^ t =b, i=l,2. 



Me- g (b)). Xl=X2= ff -i( e ) 7 fc = Ai + A 2 + Ai 2 . 



F d (xi,x 2 ) 

F c (xi,x 2 ) = e k[g(x) - g(b)] , x 1= x 2 = x, g-^e) <x<b 

and 

F ac (xi X 2 ) = e Xl ^( Xl ^~9{b)] + \2[g{x2)-g(b)] + \i2 m&x[g(x 1 )-g(b),g(x 2 )~g(b)] 

g~ 1 (e) < Xi < b; Xi > 0; i = 1,2; A12 > 0. Few members of ( 13.111 ) are listed in 
table [T] 
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Table 1. 



No 


x * y 


e 




Distribution 


1 


x + y 





g(x) = x, 

a; 6 (— oo, b), 6 < oo 


Bivariate Type 3 extreme value distribution 

F(x\, X2) = exp[Ai (xi — b) + — b) + A12 max(xi — 6, x-2 — b)] —00 < Xj < 6; Aj > 0, 
i = 1,2; A12 > 0. 


2 


xy 


1 


g(x) = logs, 
i e (0,6), 6 < oo 


Bivariate power function distribution 

F(xi,x 2 ) = ( 2 s L ) Al ( £ ^) A2 . e Al2max [ log (Tr)'(Tr)' 
< Xi < b; Xi > 0, i = 1, 2; A12 > 0. 


3 


x+y + xy 





g(x) = log(x + 1), 
x e (-1, 6), 6 < oo 


Bivariate power function distribution 

F(X 1; X 2 ) = (^)A l( ^ ) A 2e A 12 max[lo g (^),log(^ + i)] 

-1 < Xi < b; A, > 0, i = 1, 2; A i2 > 0. 


4 


■\A 2 + y 2 





g(x) = -x 2 , 
x g (— oo, 0) 


Bivariate reflected Weibull distribution 

F(xi, x 2 ) = cxp[— Aix 2 — A2X 2 , + A12 max(-ij, — x|)] 
-00 < Xi < 0; A, > 0, i = 1, 2; A12 > 0. 



Asha G. and Rejeesh C. John/Probability models characterized by 



9 



References 

ACZEL, J. (1966), Lectures on functional equations and their applications, Academic 
Press, New York. 

Andersen, P. K., Borgan, O., Gill, R. D. and Keiding, N. (1993), Statistical 
methods based on counting processes, Springer- Verlag, New York. 

Asha, G. and Rejeesh, C. J. (2007), Models characterized by the reversed lack of 
memory property, Calcutta Statistical Association Bulletin, Vol. 59, Nos. 233-234, 
1-14. 

Castillo, E. and Hadi, A. S. (1995), 'Modelling life time data with applications to 
fatigue models', Journal of American Statistical Association 90(431), 1041-1054. 

GUPTA, R. C. AND Hann, W. (2001), Analysing survival data by proportional re- 
versed hazard model', International Journal of Reliability and Applications, 2, 203- 
216. 

Johnson, N. L., Kotz, S. and Balakrishnan, N. (1995), Continuous univariate 

distributions, Vol. 2, John Wiley and Sons, New York. 
LAI, C. D. AND XlE, M. (2005), Stochastic aging and dependence for reliability, 

Springer- Verlag, New York. 
MULIERE, P. AND Prakasa Rao, B. L. S. (2003), 'Characterization of probability 

distributions under binary associative operation', Sankhya Series A, 65, 799-806. 
MULIERE, P. AND SCARCINI, M. (1987), 'Characterization of a Marshall-Olkin type 

class of distributions', Annals of Institute of Mathematics, 39, 429^-41. 
Prakasa Rao, B. L. S. (2004), 'On a bivariate lack of memory property under 

binary associative operation', Communications in Statistics-Theory and Methods, 

33(12), 3103-3114. 

Roy, D. (2002), 'On bivariate lack of memory property and a new definition', Annals 
of Institute of Statistical Mathematics, 54, 404^1-10. 



imsart-ejs ver. 2008/01/24 file: ejs_2 8_30 7.tex date: October 6, 2008 



